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GLOBAL SOLUTIONS TO THE OLDROYD-B MODEL WITH A CLASS OE LARGE 

INITIAL DATA 

DAOYUAN FANG, RUIZHAO ZI 


Abstract. Consider a global wellposed problem for the incompressible Oldroyd-B model. It is shown 
that this set of equations admits a unique global solution provided the initial horizontal velocity u^, 
the product WMg of the coupling parameter a) and initial the vertical velocity and initial symmetric 

. d 

tensor of constrains To are sufficient small in the scaling invariant Besov space Sj, x fijpd > 2. In 
particular, the result implies the global well-posedness of Oldroyd-B model with large initial vertical 
velocity u^. 


1. Introduction 


In this paper, we focus on a prototypical model for viscoelastic fluids, the so called Oldroyd-B model, 
which was originally introduced by J.G. Oldroyd ESI . This type of fluids is described by the follow¬ 
ing set of equations 


( 1 . 1 ) 


Re {ut + {u ■ V)m) - (1 - aj)Au + VII = divr, 
■ We(T( -I- (u ■ V)t -I- ga{T, Vm)) -I- t = 2u)D{u), 
divM = 0, 


where u is the velocity of the fluid, t is the symmetric tensor of constrains, and 11 is the pressure which 
can be regarded as the Lagrange multiplier for the divergence free condition. Re and We denote the 
Reynolds and Weissenberg number of the fluid, respectively, o) € (0,1) is the coupling parameter. 
Moreover, the quadratic term ga is given by ga{T, Vn) tW{u) - W{u)T-a {D{u)t + tD{u)) for some 
a e [-1,1], where D{u) = 1/2(Vm -i- (Vm)^) is the deformation tensor and W{u) = 1/2(Vm - (Vn)^) is 
the vorticity tensor. 

A brief derivation of system (11.11 ) is as follows. In fact, the dynamics of homogeneous, isothermal, 
and incompressible fluid is described by 


J dfU + {u ■ V)m = divT, 

^ ^ I divM = 0, 

where n, T are the velocity and stress tensor, respectively. Moreover, T can be decomposed into 
-p Id -I- T, where p denotes the pressure of the fluid and r the tangential part of the stress tensor, 
respectively. For the Oldroyd-B model, r is given by the relation 


(1.3) 


T + A- 


Dt 


= 2t] D{u) + p 


DaD{u) \ 
Dt j’ 


< 7 ^ 

where denotes the “objective derivative” 

= d,T + (u ■ V)t -b gair, Vm). 

The parameters A > p > 0 denote the relaxation and retardation time, respectively, and rj is the 
viscosity of the fluid. 
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The tangential part of the stress tensor t can be decomposed as r = where - 2ijjD{u) 

corresponds to the Newtonian part and to the purely elastic part. Denoting = r, the above 
equations can be rewritten as 


(1.4) 


dtU + {u ■ V)m — rj{\ — co)Au + Vp 
■ divM 

A{dtT + {u ■ V)t + ga(T, Vm)) + T 


divT, 

0 , 

2ria)D{u), 


with the coupling parameter (u:=l-^€(0,1). Using dimensionless variables, Oldroyd-B model 
can be thus described by the system (11.11) . For more details of the modeling, c.f., OTlfT^ for example. 

The theory of Oldroyd-B fluids recently gained quite some attention. According to the coupling 
strength between the velocity u and the symmetric tensor of the constrains r, the results on Oldroyd-B 
model fall into two categories. 

(I) Weak coupling case: small coupling parameter co. The study of the incompressible Oldroyd- 
B model was started by a pioneering paper given by Guillope and Saut lITSl . They proved that (i) 
system (11.11) admits a unique local strong solution in suitable Sobolev spaces H%Q.) for bounded 
domains Q c ; and (ii) this solution is global provided the data as well as the coupling parameter 
oj are sufficiently small. For extensions to this results to the L^-setting, see the work of Femandez- 
Cara, Guillen and Ortega lfT3l . The situation of exterior domains was considered first by Hieber, 
Naito and Shibata HU, where the existence of a unique global strong solution defined in certain 
function spaces was proved provided the initial data and the coupling parameter m are small enough. 
Another remarkable way to study the Oldroyd-B model (11.11) is constructing solutions in scaling 
invariant spaces. This approach goes back to the seminal work by Fujita and Kato lfT4l for the 
classical incompressible Navier-Stokes equations. While for the Oldroyd-B model, strictly speaking, 
the system (11.11) does not have any scaling invariance. Fortunately, if we neglect the coupling term 
divT and the damping term t, it is found that (11.11) is invariant under the transformation 


I (uo, To) {£uo{€x), to{€x)), 

\{u{t, x), T{t, x), n(t, x)) {iu{ix, i^t), T{ix, i^t), (^Il{ix, i^t)), 

for any C > 0. Chemin and Masmodi 0 first study the Oldroyd-B model (11.11) alone this line. For 

small coupling parameter to, they showed the global well-posedness of system (11.11) with initial data 

/ d_]\d / d \dxd 

{uq, To) belonging to the critical space Ij x IJ for any p e [1, cxj). 

(II) Strong coupling case: arbitrary coupling parameter co € (0,1). Oldroyd-B model with strong 
coupling parameter describes the behavior of the viscoelastic fluids more appropriately. The first 
result on the system (11.11) with strong coupling parameter was obtained by Molinet and Talhouk ll27l 
for bounded domains. They removed the smallness restriction on the coupling parameter co in lfT5]| . 
For the situation of exterior domains, Hieber and the authors recently ifTTIl improved the main result 
given in ifT^ to the situation of non-small coupling parameter. Constructing solutions in scaling 
invariant spaces with strong coupling is more difficult. Some efforts were made by Chen and Miao 
in fSl, where they established global solutions to the incompressible Oldroyd-B model with small 
initial data in oo’ ^ I ■ Later on, this result was generalized by the second author If34l with data in 

//^nB|j,-|<5<|-l. It was not until very recently that Zhang and the authors |[3^ obtained the 
scaling invariant solution with strong coupling in the critical LP framework. For the weak solutions 
of incompressible Oldroyd-B fluids, see the work of Lions and Masmoudi ll2^ for the case cr = 0. 
The general case a 0 is still open up to now. As for the blow-up criterions of the incompressible 
Oldroyd-B model, there are works IT) HH |22l. Besides, we would like to mention that Constantin 
and Kliegl f9l proved the global regularity of solutions in two dimensional case for the Oldroyd-B 
fluids with diffusive stress. An approach based on Lagrangian particle dynamics can be found in 
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On the other hand, for the classical Navier-Stokes equations, taking advantage of the the algebraical 
structure of the incompressible condition divn = 0, some recent progresses have been made such as 
ll^ 1^ 1511^. Following these ideas, Zhang 071 proved the existence and uniqueness of solutions to 
system related to the incompressible viscoelastic fluids with a class of large initial data. 

Motivated by O^OTl . the aim of this paper is to investigate the global well-posedness of system 
(11.11) in critical Besov spaces with the vertical component of the initial velocity field being large. Ba¬ 
sically it is the algebraical structure of system (11.11) that enables us to consider this type of initial data. 
More precisely, using the divergence free condition divn = 0, the vertical component of the convec¬ 
tion term u ■ can be rewritten as u^dhuf^ - u^diwhu^. As a result, the vertical component u'^ of the 
velocity satisfies a linear equation with coefficients depending on the horizontal components and 
T. This is the reason why we do not believe we have to impose the smallness condition on the initial 
vertical velocity u^. While the equation for the horizontal velocity contains bilinear terms in the 
horizontal components and also terms describing the interactions between the horizontal components 
and the vertical one. This necessitates the smallness restriction on the initial horizontal velocity Mq. 
Besides, since the equation of the symmetric tensor of constrains r is nonlinear, considering the inter¬ 
actions between the velocity u (with large vertical component) and r, imposing a smallness condition 
on To is necessary. In addition, we would like to emphasize that both the equation of vertical velocity 
and horizontal velocity are linearly coupled with t, and conversely, the equation of the symmetric 
tensor of constrains t is linearly coupled with u as well. In particular, the coupling parameter a> is 
involved. In order to decouple u and r linearly, one way to do is simply to apply the estimates for 
heat equation and transport equation with damping to the velocity u and the symmetric tensor of 
constrains r respectively, and to substitute the estimate for t to that for u. Nevertheless, in doing this, 
the coupling parameter oj must be small enough to close the estimate of u. This less interesting weak 
coupling case is exactly what we want to avoid considering. 

Although the strong coupling case with initial data in type critical Besov space was investigated 
by us in ll33l . large initial vertical velocity is not allowed in |[33l . A natural question is that can we 
solve system (11.11) with strong coupling parameter m and large initial vertical velocity u^l Before 
answering it, let us give some heuristic observations. Given large initial vertical velocity, it has been 
pointed out above that the horizontal velocity and the symmetric tensor of constrains r should 
be small. In the circumstances, the linear term 2a)D{u) should also be small otherwise the equation 
(fO o of T can not hold. Since we can not expect the vertical velocity to be small, it is reasonable 
to construction solutions to system (11.11) with the product of the coupling parameter m and the 
vertical velocity small. That’s why we need the initial product mu^ to be small in this paper. 

Let us now explain the main ingredients of the proof. Noting that the algebraical structure of system 
(11.11) plays a key role in our result, the method used in Il33l which relied heavily on the Green matrix 
of the linearized system of (11.11) and was carried our by Lagrangian approach fails to work here. It is 
because the divergence free condition does not hold anymore in Lagrangian coordinate. In this paper, 
we exploit energy approach in the critical framework. As the damping r in the equation (ll.lb o may 
mislead us to the weak coupling case, so a crucial step would be to find a structural mechanism which 
indicates a parabolic smoothing effect on u in the absence of the damping effect on r. To the best 
of our knowledge, in framework, researchers ETl [TTl [8l |34l were inclined to use the cancelation 
relation 


(1.6) 


(divTln) -I- (D(m)|t) = 0 


to eliminate the effect of the coupling parameter oj. The price they had to pay was that the solution 
space served for u and r would be of the same regularity, which violated the scaling. So our strategy 
is not to estimate {u, r) as a whole with an eye to (11.61) . but to find a new quantity which obeys the 
scaling and can be used as a substitute for r. An ideal candidate is Pdivr because on the one hand, 
PdivT is the exact quantity linearly coupled with u and obeying the scaling. On the other hand, it has 


4 


DAOYUAN FANG, RUIZHAO ZI 


been shown in |[^ that the Green matrix G{t, x) of the following auxiliary system of {u, Pdivr) 

Ren, - (1 - a>)AM - Pdivr = -ReP(M • V)m, 

We Pdivr, + Pdivr - cuAn = -We (Pdiv(M • V)r + PdivgQ.(r, Vm)) , 

(1-7) 

m(0) = uo, 

Pdivr(O) = Pdivro- 


behaves just like the heat kernel in the low frequency case. In this paper, we will pay much attention 
to the system (11.71) in the spirit of Danchin IfTOll . More precisely, we first localize (11.71) in frequencies 
according to the Littlewood-Paley decomposition. Then a variant energy method will be used to 
estimate each dyadic block (AqU, A^Pdivr). Owing to the fact that the linear operator associated with 
(11.71) can not be diagonalized in a basis independent of coercive estimates can not be achieved 
by standard energy argument. Since the parabolic-hyperbolic system (11.71) behaves differently in high 




\\u, 


q\\]2 + 2 


_ II 


2\\Uq\\^ 




(l-^)We p^. ||2 _ 


+ 2(M^|^Pdivr^) -I- 2(M^|^APdivr^) 


IIg 



to be chosen later) 

for 

q < qi. 

for 

q > qo, 

for 

q\ < q < qo- 


Some words about the dehnition of Yq are in order. 

• For q < q\, the parabolic smoothing effect of Pdivr is obtained with the aid of the linear 
coupling term -Pdivr in (ll.7l) i. that’s why we introduce the cross term (u^l^APdivr^). 

• For q > qo, compared with the linear coupling term -Pdivr in (11.71) ] . the damping term Pdivr 
in (0)2 is predominate, and thus there is a damping effect on Pdivr. 

• For q\ < q < qo, a. new quantity A“^Pdivr is introduced to hll up the gap between qi and qo. 
Indeed, in this case, we do not need to care about the scaling thanks to Bernstein’s inequalities. 
On this basis, the cancelation relation (11.61) motivates us to choose A“'Pdivr to replace Pdivr. 
In doing so, u and A“'Pdivr are linearly decoupled, and A“^ Pdivr also possesses a damping 
effect because of the damping term. Therefore, in this paper we combine the case q > qo 
with qi < q < qo, and regard them as high frequency. The left part ^ is regarded as low 
frequency. 

We would like to remark that similar construction of Yq was also used by the second author in ll^ 
for the compressible Oldroyd-B model. However, for qi < q < qo, the quantity in OSl corresponding 
to Yq is 




q\\l2 


+ 


We 

(uRe 


l|A- 


■ipdivr,||2, 


2y3(l - m)We 
mRe 


(PwglPdivr^) 


with some positive constant /3. Since the cross term (Pu^lPdivr^) is involved in Yq, A“^Pdivr in |[35l 
exhibits the parabolic smoothing effect stemming from the linear coupling term -Pdivr. Hence the 
case q\ < q < qoi^ regarded as low frequency in IMl . The point is, though, that Yq can not be used 
in this paper to replace Yq because the parabolic smoothing effect for A“^ Pdivr and hence for u is too 
weak (depending on a») to deal with the linear coupling term 2a)D{u) in (II.H o. 

Obviously, the above analysis do not take into account the anisotropy. In order to distinguish the 
velocity u from different directions, we must bound {u!\ (Pdivr)^‘) and (m^, (Pdivr)'^) separately (please 
refer to Section [3]for more details). 

We shall obtain the existence and uniqueness of a solution {u, r) to (11.11) in the following space. 


Definition 1.1. For T > 0, and 5 € R, let us denote 

■.= n 4(R^y)/ X (Ct(4]) n . 
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'We use the notation if T = oo, changing [0, T] into [0, oo) in the definition above. 


Our main result is stated as follows: 


d , \d 

f I 


d \dxd 

>2 


Theorem 1.1, Let d > 2. Assume that {uo,to) e i " j ^ {^2 ij whh divuo = 0. There exists a 
positive constant M, depending on d, Re and We, but independent of co,such that if 

\ 


( 1 . 8 ) 


Co 


- -^ IImqII . d , + oj\\uq\\ 1 + IItoII , exp 

O V ^ 2,1 ^ 2,1 ) 




(1-m) 


12 "“oil .4-1 


< 1 , 


then system (O) admits a unique global solution {u, r) in £2 with 

llu'^IL . 4 _, +||(PdiVT)^|L . 4 _, +||u^|| , . 4 ,, + ||(PdiVT)^|| ,d d 

) ir(B 2 i ) ^ 1(^21 ) C (®2 1 ) 


(1.9) 

and 


< 


Co 


(1 - m)3 


1 + ||moII. 4 _i + ||(PdivT)ol| 4 _i I, 

^2A ^2,1 


( 1 . 10 ) 

and 


\\u'^\\_ . 4 _„ + ||(PdiVT)'*|L . 4 _. +\\u% , . 4.1 + ||(PdiVT)'’|| d,,d 

LTiBf ) LfBf ) L\(Bf ^ ) 


Co 




IImqII . d , + mllwoll. d , + llroll. d^ \ exp 




( 1 -W)I 2 " ‘"'si-' 


( 1 . 11 ) 


< 


LT(Bf) 

Co 


.4 +lml , .4 

r 2 \ 2 Lr 2 


L](BIP 


(l-m )6 


ll“oll. 4-1 + ^ll'^oll ■ 4-1 ll’^oll. 4-1 
^21 ^21 s|, y 


exp 




(1-m) 




The rest part of this paper is organized as follows. In Section 2, we introduce the tools ( the 
Littlewood-Paley decomposition and paradifferertial calculus) and give some nonlinear estimates in 
Besov spaces. Section 3 is devoted to the global a priori estimates of system (11.11) . The proof of 
Theorem II.Il ls given in Section 4. 


Notation. 

( 1 ) ¥0Vf€S’,fq-kj. 

(2) Forx = (xi,--- ,Xd), we denote x/j := , x^-i); For vector field m(x) = (m^(x), •• • ,m^(x)), 

we denote u^{x) := (u^{x), • • • , u^^^{x)), dhu{x) {diu{x), ■ ■ ■ , dd-iu{x)y, and div/,M^(x) 
d\u^{x) + • • • + dd-iu‘‘~^(x). 

(3) For a,b € L^, {a\b) denotes the inner product of a and b. 

(4) Throughout the paper, C denotes various “harmless” positive constants, and we sometimes 
use the notation A < B as an equivalent to A < CB. The notation A ~ B means that A < B 
and B < A. 


2. The Functional Tool Box 

The results of the present paper rely on the use of a dyadic partition of unity with respect to the Fourier 
variables, the so-called the Littlewood-Paley decomposition. Let us briefly explain how it may be built 
in the case x € which the readers may see more details in [[T] |3l. Let ix, ip) be a couple of C°° 
functions satisfying 

Supp.^ c (1^1 < ^}, Suppi^o c {^ < 1^1 < ^), 


(2.1) 
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and 

x{^)+ Y,^(2-‘^^) = I, 

q>0 

for ^^0. 

qeX 

Set ipq{^) = hq - T'~^{ipq), and h = T~^(x)- The dyadic blocks and the low-frequency cutoff 

operators are defined for all g e Z by 

hqU — ^(2“^D)m = 

=;^(2-^D)n- 

Then 

(2.2) u = ^ A^m, 

^€Z 

holds for tempered distributions modulo polynomials. As working modulo polynomials is not appro¬ 
priate for nonlinear problems, we shall restrict our attention to the set of tempered distributions u 
such that 

lim ||S^m||l»=0 . 

q-^-oo 

Note that (12.21) holds true whenever u is in S'j^ and that one may write 



p<q-l 


Besides, we would like to mention that the Littlewood-Paley decomposition has a nice property of 
quasi-orthogonality: 

(2.3) ApA^M = 0 if Ip - ^1 > 2 and Ap{Sq-iuAqu) = 0 if |p - ^| > 5. 

One can now give the definition of homogeneous Besov spaces. 

Definition 2.1. For s eM., {p, r) € [ 1 , oo]^, and u € .S'(R^), we set 

\\uUs^^^ = \\VnAqU\\Lp\l,. 

We then define the space Bp ^. := {u e IImIIb* ^ < oo). 

Since homogeneous Besov spaces fail to have nice inclusion properties, it is wise to define hybrid 
Besov spaces where the growth conditions satisfied by the dyadic blocks are different for low and high 
frequencies. In fact, hybrid Besov spaces played a crucial role for proving global well-posedness of 
compressible barotropic Navier-Stokes equations in critical spaces ifTOll . Let us now define the hybrid 
Besov spaces that we need. Here our notations are somehow different from those in ifTOl . 

Definition 2.2. Let 5 , t € R, and u e 5'(R^). For some fixed qo € Z, we set 

\M\b-\ ■= 2] 2‘?^llA,n|lL2 + 2] 2^1|A,n||p2. 
q<qo q>qo 

We then define the space := {u e .Sp(R^), llnllgw < oo). 

Remark 2.1. For all s,t eM.,qQ e Z, and u e .S'(R^), setting 

llnllgv := 2] 2^^||A,n||p2 + 2] 2‘^’\\Aqu\\p2, 

q<qo q>qo 

then it is easy to verify that ||m||^.w w ||M||g.w . 

The following lemma describes the way derivatives act on spectrally localized functions. 


I hq(y)u{x - y)dy, 
I hq(y)u{x - y)dy. 
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Lemma 2.1 (Bernstein’s inequalities). Let A: € N and 0 < r < R. There exists a constant C depending 
on r,R and d such that for all (a, b) € [1, oo]^, we have for all A > 0 and multi-index a 

. //Supp/c then sup„=, ||5 "/IIl^ < A^^‘i^'a-Tt)\\f\\^.. 

. //Supp/cC(0,^r,^/?), then C-^-^A^\\f\y < sup|„|=, ||5 “/IIl“ < 

In fact, the following special case of Lemma IZT] will be used very frequently in this paper. 

Corollary 2.1. Let A V-A, then for any f € S', there holds 

(2.4) \2%\f\\L^ < IIAA,/||l2 < ^-in'AJWiy. 

Owing to the property (12.11) of Supp(/), (12.41) can be easily obtained by virtue of Plancherel’s Theorem. 
Let us now state some classical properties for the Besov spaces. 


Proposition 2.1. For all s, sy, S 2 1 < p, P\, P 2 ,r,r\,r 2 < oa, the following properties hold true: 

• If P\ < P 2 and ri < r 2 , then ^ Bp^ ^l . 

• If Si S 2 and 9 € (0,1), 

• For any smooth homogeneous of degree m e "Z function F on R'^\{0), the operator F(D) maps 


Next we recall a few nonlinear estimates in Besov spaces which may be obtained by means of 
paradifferential calculus. Firstly introduced by J. M. Bony in 0, the paraproduct between / and g is 
defined by 

ffg-Y,^,-ifA,g, 

qeX 


and the remainder is given by 

B{f,g) = ^ A^/Agg 
q>-l 


with 


Ayg = (Ag_i + Ag + Ag+i)g. 

We have the following so-called Bony’s decomposition: 

(2.5) fg = tfg + tj + R{f, g) = tfg + t'J, 

where t'^f := tgf + R{f,g). The paraproduct t and the remainder R operators satisfy the following 
continuous properties. 


Proposition 2.2. For all s cr > Q, and\ < p,p\,p 2 ,r,ri,r 2 < oa, the paraproduct t isabilinear, 
continuous operator from L°° xBp^toBp^ and from B'^^^ X Bp^.^ to B^~f with 4 = min{l, ^ 

The remainder R is bilinear continuous from Bp\ri X ®p 2 ,f 2 Bp\r^^ with si-t- S 2 > 0 , 2 = ^ 

and 4 = — + J- < p 

r ri r 2 - 

In view of (12.51) . Proposition 12.21 and Bernstein’s inequalities, one easily deduces the following 
product estimates: 

Corollary 2.2. Let p € [l,oo]. If si,S 2 < ^ and -i- ^2 > dmax{0, | - 1), then there holds 

(2.6) ||mv|| rf < C||M||g.-i ||v||g^2 . 

B . P.i P.l 

In the following, we shall give a commutator estimate, which will be used to deal with the convec¬ 
tion terms. 


Lemma 2.2 ( ll35l ). Let A := V-A, then 

(2.7) ||[A-i,A,_iv V]A,n||i2 < C||VS,-iv||l»||A-'A, n||^2. 
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The divergence free condition allows us to get the following lemma, which plays a key role in this 
paper. 


Lemma 2.3 (|[37l). Let V be a divergence free vector field. We then have 


(LT) ^ C\\dmv\2\\vX2, 


and IIv^IIl” <C||v'*||^ llv^f. 






The study of non-stationary PDEs requires spaces of the type LjiX) - L^(0, T-,X) for appropriate 
Banach spaces X. In our case, we expect X to be a Besov space, so that it is natural to localize the 
equations through Littlewood-Paley decomposition. We then get estimates for each dyadic block and 
perform integration in time. But, in doing so, we obtain the bounds in spaces which are not of the 
type L^(0, T ; Bp ,.). That naturally leads to the following definition introduced by Chemin and Lemer 
in 14]. 


Definition 2.3. For p € [1, +oo], € R, and T e (0, +oo), we set 


II«IIk 




2 ^'*||A^M(t)||^P(^) 


and denote by Lj{Bpfi the subset of distributions u e 5'([0, T] x R^) with finite \\u\\iPps^^^) norm. 
When T = +oo, the index T is omitted. We further denote CjiBp f) = C([0, T]; n Lf{Bp fi. 


Remark 2.2. All the properties of continuity for the paraproduct, remainder, and product remain true 
for the Chemin-Lemer spaces. The exponent p just has to behave according to Holder’s ineauality 
for the time variable. 


Remark 2.3. The spaces Lp{Bp can be linked with the classical space Lj{Bp via the Minkowski 
inequality: 


r(Skr) 


< 


M\Li( 


if r>p. 


IImIIz 




if r < p. 


3. A Priori Estimates 


We begin this section by rewriting the system (11.11) in term of {u‘, (Pdivr)'), 1 < / < d, as follows 

fc-^Au'-^(PdivTy = /', 

(Pdivr)' + ^(Pdivr)' - - h\ 

where f -(P(m • V)m)', and H := - (Pdiv(M • V)t + PdivgQ,(T, Xu))'. Eor simplicity, let us denote 
cr := Pdivr, then system (13.11) can be localized as 


(3.2) 


d,Uq +Sq-lU- Xu\ - ^Au'q - - f^, 

dtO-'q +Sq-lU- Xo-\ + - ^Au^ = h‘^. 


where f^ := Sq-iu ■ + /^, and h'^ := Sq-\u ■ Xcr'^ + h'^. 

Denote 


a'Xo 

A\t) 

a\0) 

Bit) 


IKIL 


L-TiBf 


- 1 ) 




+ inL ... +l|u'' 


,4-1 +M- .4-1 +IK 


Ll(Blf) 


.40 + 11 ^^ 


LTiBf ) 


LTiBf ) 


LfBf ) 


LfBf ^ ) 

+ lki , .d u-i, 

L\(Bf ^ ) 


ll^oll. 4-1 + IKoll . 4-1 ’ A'^(O) ||mq|| . . + ||crQ|| d^ , 

®2,1 ®2,1 ^21 ^21 


IL .d + 
lT(bIp 


l){bIp 


. , B(0 ):-||ro||... 


We will establish the following proposition. 
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Proposition 3.1. Let {u, t) be a solution to system (11.11) . There exist constants Ci, C 2 and C 3 depend¬ 
ing on d, Re, and We but independent ofco, such that the following inequalities hold 


(3.3) 


A\t) < 


‘a\0) + 


(1 - 09)3 


f 

Jo 


+ 


. S + 1 


dA\\T\\,ddt' 


^*".9 +||«^||^ ||«^||^ ||«^||^ ||«''||^ 

r2 r2 + ‘ aT + I 


B, 


2,1 


2,1 


B. 


2,1 


B. 


2,1 


^ 2 “l / 


dt' 


f 


,hii 2 




d 


+ 1 " ". 4+1 




2,1 


A‘^(t) < 


'a^(O) + 


(1 - 09)3 


d^, + 




"^ 2,1 


.d\\l 


iJWrW.ddt' 




+ 


f 


(3.4) 
and 

(3.5) B{t) < C 3 1 B{0) + CO 




hu +ii„ft||5 


■4+ 

Bl 


h\\2 


d _ 

Bl 


dll 


d _ 


d 

Hi 


dt' 


■ ir-t . 

Bl 


d\l 


Bl 


Hi 9 


d-idt 1 , 


rfj-l + 


L}(Hi j 


LjiHiJ 



. 4 . 1 + (II«'’II.4.i+II“''II.4.i)IMI.4^^1- 


- ’T + 

Bl 


Hi 


Proof. For \ <i <d, let us denote 


Fi 



We , 

\\u%+2 



.We , 


We 




+ 2(u\\—(T'^) + 2{u\\—Aa\) 


, for q< qi. 


Y’ ■= 


2\\uX2 + 


(1 - 09)We , 




,,(l-09)We 
09Re 


, for q > qo. 


y' — 

^ • 






, for qi <q< qo. 


where ^0 and qi are to be determined later. 


Estimates of F^, q eZ. 

Step 1: low frequencies.-q < qi. 

Taking the Lf inner product of (13.21) ; and (l3.2l) o with u‘^ and cr'^ respectively, using the divergence free 
condition divu = 0 , we obtain 

(3-6) 

and 


(3.7) = ^KK)- 

To cancel the “bad” cross term (u^lcrp in (13.61) . we need the following equality for the Lf inner product 
(m^Ict^) of and cr^, 


CO 


1 


(3.8) -«k;) + —IIAu'11^2 - —\K\\r2 - 


1 


dt 


We 


Re 




-^{Au‘^\(r'g) + :^{u‘^\cr‘^) - (/^Icrp + 


■'qi-'-qJ 
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which is nothing but the sum of ^ (I3.2l) i and (^ (I3.2 I) tIm 0. Multiplying (13.81 ) by adding the 
resulting equality to (13.61) . one deduces that 


1 d 

2 dt 


We 




1 

Re 


^llL 2 


We 


\\(r% 


(3.9) 


1 -m ■ We . . . We/ • • • • x 


Re 


Multiplying (13.71 ) by 2 , adding the resulting equality to (13.91 ). then the negative term ~ 

in (13.91 ) is absorbed, and we have 


2 


1 d 

2 dt 


lKlli2+2||^cr;||^,+2(n;i^crp 


,We ; 2 ;.We 




1 + m ,■ We 


Re 


(3.10) 


• • (WeY We / • • • • x 

+ 2 ( 1 ^) w,K-> +((/.;ity)+ (kk^) 


In order to exhibit the smoothing effect of cr', 1 < / < d in the low frequency case, we shall make use 
of the linear coupling term -^cr' in the first equation of (13.11) . To do this, taking the inner product 
of (I3.2l) i with Ao"^, and taking the inner product of (l3.2l) o with Am^, then adding them together, we 
are led to 


|("ilAcr',) + T||A,ril|J. - ,^||A„i"2 


i^(An;|A(r') + i^(An>') 


(3.11) = (/'|A(r') + (^'|Anp-(S,-in-Vn^|A(rp-(,S,_in-V(r'|Anp. 

Multiplying (13.111) by ^, adding the resulting equality to (13.101) yields 

1 d 

2 dt 


■ o We • • We • We 

+ 2|l^cr;||2, + 2(n;|—crp + 2(n;|—Acr^) 


1 


--^(Anil^Atri) + f J-- 


OJ 


.i l|2 


1 ..We 


i ii2 


1 ..We 


111^2 


Re 


We 


Re / 




Re 


“^'Re 

• • /We\^ We / ■ • • • X 

W + 2 ( 1 ^) + 157 (ifiK) + %V,)) 


Re 


(3.12) 


+ 


We 


((/glAcrp + (hqlAUg) - {Sq-lU ■ VUglAo-'g) - {Sq-lU ■ Vo-'qlAu'g)) , 


Let us now estimate the cross terms contained in T^ for q < q\. Indeed, using Cauchy-Schwarz 
inequality and (12.41) . one easily deduces that 


2 |(m' l^(z')| < ^-\\u% + || 2 „ 


and 


.We 


Re 


2 |(n:,|—A ct')| < | - 2 ^‘ | | \\u%, + ||—(t; 


.We 


Re 


411^2 


1 / . 0 We . 0 
^ 4 + Il^‘-;ili 2 


provided 

(3.13) 

Consequently, 




We 


•• ..o 


We 


2 '?* < —. 

“ 16 


- - We, 


;• 1,2 


3 ..We 


i 
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II 


and hence 
(3.14) 


+ if q<qu 


with <71 satisfying (13.131) . Thanks to (12.41) . the remaining two cross terms on the left hand side of 
(13.121 ) can be bounded as follows: 


1 

We 


1 + m 
Re 


We 


^Re + We 8 _,„^ We 
- ^ wp ^2‘i'\\Ku \\i2 —cr 

WeRe 3 ^ Re ^ z ,2 

8 , Re + We / 1 , 1 


We 


(3.15) 
provided 

(3.16) 
Moreover, 




We , 


L 2 


2 ^‘ < 


3 VWeRe 
32 Re + We’ 


1 - m 
Re 


We 

(An^l—Acrp 


^ k2"‘ IIAn; 


(3.17) 
provided 

(3.18) 

Finally, using (12.41) again yields 


2 \3 
1 1 


Re 


“ 8 Re2 We 




+ 


We. 

—Acr 
Re ^ 


L 2 


1 1 


We 

—Acr 

Re 


L 2 


2 ^‘<—, and 2 ‘?‘<-J—. 
16 8 VWe 


(3.19) 




||2 < 1 12'?' 1 ||An'||" < A^||An'||" , 


1 


ii 2 


1 1 


Re 




4 Re' 


Re" 

provided (13.131 ) holds. Now collecting (13.131 ). (13.161) and (13.181 ). we take 

3 VWeRe 'l 




(3.20) 


q\ ■= 


log2 


32 Re + We 


From (I3.15I) - (I3.19I ). and using the left hand side of (12.41) . we find that, for q < qi with q\ satisfying 
(13.201) . there holds 


1 




a> 


.i i |2 


1 ..We 


j i |2 


1 ..We 




1-m^^ We / 1 ,-,We ,. 

—(Am„ — Ao-^) + —--— (Am„ — o-„ 


Re 


Re 


3 1 


“ 8 Re2 We 


1 1 


(3.21) 


3 /3 


8 \4 


1 1 


We 


We 

—Acr 
Re ^ 


^ min-j^,-^!>(||M(||22 + 


We Re 




Re / 

^ +- — 

l 2 4 m 

We 


We 


Re 


L 2 


-cr. 


) 

2 

L 2 


To bound the right hand side of (13.121) . firstly integrating by parts and using the divergence free 
condition divn = 0 yields 

- {Sq-iu ■ Vu' lAcr') - {Sq-i ■ Vo-‘\Au‘) ^ (Sq-id jU^\dku‘djO-‘ + dkcr' dju'). 


(3.22) 
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Accordingly, using Bernstein’s inequality, (13.201) and Lemma 1231 we are led to 
We 


— |(S^-1M • Vm^IA cr^) + (S^_i • Vo-^lA m^)| 
^ We ; 1 

which, together with (13.141) . implies that 
(3.23) R.H.S.of (ITT2]) < C 


h: 


We- 

1141b 


" .rf+l 


,hu2 


We 

\K\\U + 


1 




Y' 

r 


Substituting (13.231) into (13.121) , we infer from (13.141) and (13.211) that there exists a constant C depend¬ 
ing on d. Re, We, but independent of a», such that if ^ < ^i, there holds 


(3.24) 


-(f') +22?(f') <CF' + + 


h\\ 2 






Step 2: high frequencies.-q > q\. 

Part (i). q> qo. 

First of all, we would like to point out that, in the high frequency case, the cross term -^ (An^ |cr^) 
in (13.71) is “bad”. We will cancel it with the aid of (13.81) once more. To this end, multiplying (13.71) and 

(13.81) by and respectively, then adding the resulting equalities together yields 


1 d 
IJt 

+- 


(1 - m)We 

- Cn 

wRe 


L 2 


2 m' 


(l-m)We 

cr‘ 


1 


(3.25) 


(1 - m)We 

(1 - m)We \^ 

mRe I 


(1 - m)We ,• 

- O'r, 

mRe ^ 


a>Re 


1 


1 - 0 ). 


Re 


-IIAm'II 


% 


(4I<)- 


l 2 We 
(1 - m)We 


(1 - m)We ; 

ZiRe 


mRe 


((4K) + (4K))- 


To absorb the negative term -^^||Am ^||^2 (13.251) . we multiply (13.61) by 2, and add the resulting 

equality to (13.251) to get 

2 


1 d 




(1 - m)We 
mRe 


L 2 


■ 2 |n; 


(1 - m)We 
mRe ^4' 


\- OJ. 


Re 


■||A<II 


^114 


1 


+- 


(1 - m)We 


(1 - m)We 
mRe 


(3.26) - 2(4|n;) + '^^ 


It is easy to verify that 

(3.27) 2||n' 11^2 + 


mRe 


(1 - m)We 
mRe 




(1 - m)We 


(1 - m)We 
mRe 


mRe 


((4kp + • 


- 2 m 


L 2 


(1 - a;)We 
mRe 


11411^2 + 


(1 - m)We 
cuRe 


2 

L 2 ' 


Using Cauchy-Schwarz inequality and (12.41) yields 

(1 - m)We 


2 • ■ 1 

— (cr'Ju') H-I u' 

Re ^ ^ We' " 


mRe 


1 -F m 1 
1 


(1 - m)We 

-O’iv 

wRe ^ 


1 1 
< -- 


1 


2 1 - m We 


(1 - m)We 

CTn 


mRe 


-F 


l 2 (1 - m)We 


11411^2 
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< -- 


1 1 


1 


(3.28) 
provided 

(3.29) 

Now let us take 

(3.30) 


2 1 - w We 
1 


(1 - ai)We 


1 1 
< - 


2 1 - w We 


wRe 
(1 - w)We 


-cr. 


2 (4 

+ - 2 -® 


wRe 


-cr. 


\3 

2 1 1 -6U 

[2 2 Re 




(l-ai)We'"'“^"^^ 


1111 ^ 2 , 


4 

l2-‘?o 

3 


2 ^ 1 1 -w 

(1 - w)We “ 2 Re ’ 


i.e. 


2 ?o > 


3(1 - co) 



90 


log 2 


3(1 -w) 

Then it follows from (13.281 ). (13.291) and (12.41 ) that 



+ 1 . 


1 




Re 


(l-w)We 


a>Re 

(1 - m)We 


-cr. 


L 2 


a>Re 


-cr. 


> - 


1 1 - m , 1 1 

2 2(1- m)We 


(1 - a;)We 


oiRe 


-cr. 


L 2 


> -- 


1 1 -m /3 


2 Re 


1 


\Kfr2 + - 


1 


(3.31) 


> - 




^ '2(1- m)We 

(1 - w)We 


(1 -m)We 


mRe 


-cr. 


-cr. 


wRe 

Finally, the right hand side of (13.261) can be bounded as follows, 


(3.32) 


R.H.S.of(Il26l)<C ||/'|b + 


(1 - m)We~,. 


cuRe 


L 2 


\K\\l^ + 


(1 - m)We 


mRe 


-cr. 


L 2 


Substituting (13.32b into (13.26b . using (13.27b and (13.31b . we find that there exists a constant C indepen¬ 
dent of d, Re, We, and a>, such that if ^ > qo, there holds 


(3.33) 


- 

dt V 


1 


(1 - m)We 


(nf < CT; II4IL2 + 


(1 - m)We~,. 


(uRe 


L 2 


Part (ii). qi < q < qo- 

From the definitions of q\ and q^ (see (13.20b and (I3.30b i. it is obvious that there is a gap between 
the low frequency part q < q\ and the high frequency part q > qo since qi < qo. To overcome this 
difficulty, it is necessary to resort to some new observation. Considering that ||A“^o -^||^2 is equivalent 
to ||cr^||i 2 if < q < qo, we can use ||A“ V ^||^2 instead of ||cr^||i 2 to fill up the gap. To do so, applying 
the operator A“^ to (l3.2b T. and then taking the inner product of the resulting equation with A“^cr^, 
we obtain 


i-||A-V'"2 

2dt 


?iIl 2 


-I- —||A V'11^2 ~(c<fl|cr') 

We ^ We 


(3.34) 


^ (A %\A V')-(A • Vcr')|A V'). 


- 1,^1 


- 1 ,^'' 


Multiplying (13.34b by then adding the resulting equation to (13.6b . the cross term (m'| cr') is can¬ 


celed, and we arrive at 


^4fll“oll?2-^l|A 


2 dt 

^ (faWj + 


■ cuRe 
We 


Re 


"^11^11^2 + ^I|A“V' 11^2 


cuRe 


((A-i^;|A-V^) - {A-\Sg-iu • Vir')|A-V')), 


(3.35) 
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In view of (12.41 ) and (13.201 ). it is not difficult to verify that 


1 - (0 
Re 




1 - a» /3 


Re \4 




Il 2 


cuRe 


(3.36) 




> 


2i6(Re + We)2 


Using the divergence free condition divn = 0, one easily deduces that 


(A • VcrplA Vp = ([A \ Sq-iu ■ V]o-'g)\A V^, 

then Lemma |2!2] and Lemma |23] imply 

(3.37) |(A-'(5,_in • V(r')|A-Vp| < C||VS,_in||L<»||A“< C\\u’^\\\ ^\\u\\\ JIA'V'H^,. 

% Bl 

Substituting (I3.36I) and (I3.37I) into (I3.35I) . we find that there exists a constant C depending on d, but 
independent of Re, We, and a>, such that if < q < qo, there holds 


n 3 Sf - (rf + 3^We(l - m) / .x 2 ^ • 

^ ^ dt 2 i 6 (Re + We )2 i “ 1 


We 


II4IIU + V—||A-^/j;|L2 


^ 2,1 ^ 2,1 2 


The smoothing effect of u. 

Now we are in a position to show the smoothing effect of the velocity u. In fact, in the low 
frequency part, the smoothing effect of Uq has been obtained in (13.241) . More precisely, a direct 
consequence of (13.241) and (13.141) gives 


||“^(0||^2 + 1^9(011^2 +2 ^|||m9||^1(^2) + IK^IIl/cU) 

< c( ||<( 0 )||^, + ||cr'(O)ll^, + 


(3.39) 


■ fii«''ii^,jNi^,,(iKiiu + IKIL2)H 

-'0 b}. / 


for q < q\. Nevertheless, in the high frequency part, we infer from (13.331) . (13.271) and (13.381) that 


IImUOIIlj + 


(l-m)We . 


(3.40) < C ^114(0)11^2 + 

for q > qo, and 

II“9(0 IIl2 + 


a»Re 
(1 - m)We 


1 


l 2 (1 - m)We 






(1 - m)We 


■cr„ 
mRe ^ 


L‘(U) 


<.Re 


L 2 


+ \\fq\\L](L'^) + 


(1 - m)We 


mRe 


h\ 


L‘(U) 


We 


mRe 


||a“V4o|| 


3^We(l - u) 


2i6(Re + We)2 ^ 


We 




)(U) , 


< c 


We 




(3.41) + ||»' 






We 




mRe 


for q\ < q < qo- Recalling the definitions of qi and qo, and using (12.41) again, one deduces that 

e 2 ^ 


(3.42) 


lk:,llL 2 < :r 2 ^°||A-Vl|L 2 < 


32(1 - m) V We 


Re 

^l|A"V'||i 2 , 
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and 

(3.43) \\cT\\y > ^2?‘||A-V'||^2 > 


32 

“S'- A cr' , 2 . 

28 Re + We 


By virtue of (13.421) . (13.411) implies 

3^(1 - m)We 


II“o( 0 IIl 2 + 


2 ^ 


+ 


i(')IL= 

3^(We)2(l - a)f 1 


cr„ 


< C 


34We(l - oj) 

2i6(Re + We)2''“^''^‘(^") Vm223(Re + We)2 Re 


(3.44) 


11 ^ 5 ( 0 ) 1^2 + 

/' 


mRe 


||A-V'(0)||^, + ||/'||^.(^2) 


+ 


, I 
, h\\l 


d 

Hi 


+1 . ^+1 
Hi V 


\WM + 


We 

mRe 




ellL 2 


We 

cuRe 

dt’ 




for qi < q < qo- It can be seen from (13.401) and (13.441) that we have not yet obtain the soothing effect 
of M 5 if ^ > q\. Fortunately, the damping effect of (Tq is obtained both in (13.401) and (I3.44I) . This 
enables us to get the smoothing effect of Uq by means of the equation of the velocity u. To this end, 
taking the H inner product of (I3.2l) i with integrating by parts, we have 


1 d 


i Ii 2 


+ 


1 - (U, 


II^“?IIl2 = ^(c^5l“P + HqWq)- 


Idt" ^"2-" ■ Re 
Using the left hand side of (12.41) . we easily get 

,2 


Re 


(3.45) 


II“ 5 ( 0 IIl 2 + 


1 -m /3 


Re 4 


t 2 ‘? 


5IIl;(L2) < \Wq{0)\\i2 + —l|cr^|lLl(L2) + ll/ 5 llz,l(z, 2 ). 


Multiplying (13.401 ) by 2m, and adding the resulting equation to (13.451) . we arrive at 

,2 






Re " Re 

(3.46) < C (114 (0)11^2 + ||cr;(0)||^, + 11411^.(^2) + 

for q> qo- Multiplying (|3.44D by 
(13.431) . we obtain 




2 ?IIl;(L 2 ) + ir?llLl(L 2 ) 


1 


^iil;(l 2 ))’ 

, and adding the resulting equation to (13.451) . using 


i|“«( 0 IIl 2 + .,.1^ —— lk 5 (o|| 


1 - m /3 


Re 


< 


C 


(1 - m)^ 


34ReWe(l - m) 

(ii4(®)iiu + Ik5(0)|lz,2 + ii4iiu‘(U) 


1 


4^1 iKiiu(U)-^IKIL;(l2) 


( 3 . 47 ) +11414.(^2) + fwHwKjuwK,, (ii4iil 2 + lk^lL2) Jf')’ 

-'0 Hi Hi 

for q\ < q < qo. Now collecting (13.461) . (13.471) . and (13.391 ). we conclude that 

.4_, + (1 - m)||cr'||_ 4_i +(l-m)||M'|| d^^ + ||o-'|| 


C 




l;(b 2 , ) 


LliB. 


+ 1 , 4-1 


2,1 


(1 - m )2 


o'Ut-i Ir 0 


° Jo 


+ 


f 


hill 


. IT +1 

^ 2,1 


4-1 + r 1 - 4 - 1 W 

>2 " " d 2 I 


2.1 


2.1 


(3.48) 
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with some positive constant C depending on d, Re and We, but independent of m. 

Estimates of (II 4 IL 2 + \\h\\y). 

Recalling that := S q-\u ■ + /^, and f := -(P(m • V)m)‘ = -(d,m + diA~^dm)u^dku'”, 

order to estimate Zgezit suffices to bound q-iu^dkUq\\i 2 and \\{5im 

diA~^dm)u^dku"'\\ rf_|. 

Estimates ofY,qez q-iu^dkUq\\i 2 . 

• Case 1: 1 < / < d - 1. 

Y^2^i-^^\\Sq-,u’^dku‘qh2 

qeZ 


• Case 2: i - d. 


< C\\Sq-lUh^Y,^ku%2 


< C||n||L~||n'’|| , < C 

« 2,1 


: h \\2 


+ 


Al) 


n 2 


2,1 


for 1 < / < d. 


'" 2 , 1 ; 


Y^2^i-^\\Sq-,u^dk4\\L2 < 

qeZ qeZ qeZ 

< C\\u%^\\u‘^\\ , +C\\uX«\\dmu’^\ 4 _, 

^21 ^21 

< Clln'il 4 lln^ll. . 

^ 2^1 -e?! 

Combining these two estimates, we are led to 

/ 

(3.49) 2]2(^1)||S,_i/4<IIl2 < C \\uY,, 

qeZ ®2 

Estimates o/||(d,m + dih2'^dm)u^dku'”\\ d_^. 

Using Proposition 12.21 we have 

m,n + diA-^dm)T^,dku'"\\ d, 

Bh 

< y ll(d™ + dih-^dndt^kdku'^w d__, + y Ii(da + diK-^dd)t^,dku‘^\\ d__, 

p 2 hi 

\<m,k<d-\ 2,1 l<i<d-l "2,1 

+ y ll(d™ + diA-^dm)t,dddu’"\\ d + ||(dy + diA-^dd)t,dddu‘^\\ d 

\<m<d-\ 2,1 "2,1 

/ill im ../ill 


\ 


< C|||n'1|L»||5ftn"||.4_, +||n'l|L»||5ftn^||.4_, +||n^||z.~||5dn1.4_. +||n‘'||L~l|div,n'*||.4_, 

Sza Bl Bl j 


( 


< C 


\ 


'‘" 2 , +||n^|| a lln^ll a 


r5 


hi hi 


" 2 , 1 ; 


and 


\\{6im + diA-^dJT'M.U 


< Yj ll(<^™+^i-^“^^'ii)^L"“^ll.fi + Z md + diK-^dd)t'^^^duH.d-, 


\<m,k<d-\ 


2,1 \<k<d-\ 


^ 2,1 ” 


l<m<d-l 











OLDROYD-B MODEL WITH A CLASS OF LARGE INITIAL DATA 


17 


< C ||5,m'’|| . + |Ml. 4 -JKII. 4 + ||div/y||, 


< c 


“ 2,1 “ 2,1 “ 2,1 
( \ 


<1 


^ 2,1 


2.1 


^ 2,1 


2,1 


2,1 


^ 2 , 1 / 


Therefore, according to Bony’s decomposition, there holds 
(3.50) 11(5™ + diK-^d,n)u'^dki^^\\ j,<C 






. w + - 
® 2,1 




2,1 


Ms 


for \ <i <d. 


It follows from (13.491) . (13.501) and Proposition 12. 1 I that 


< c[ 


. ,line’ll i., +\\A\ i WuH i 




'^ 2,1 


^ 2,1 


Ms 


(3.51) 


< C 




rf ||2 


V 2,1 2,1 

d 


• 4 + 1 

Ms 


h\\2 


■4-1 

Ms 


d\\2 


■ 4-1 
Ms 


Ms ! 


for \ < i < d. 


Next, we turn to estimate 2 ^*^^ 0 ||/j^||^ 2 . Noticing that 
hq = Sq-xU-Vcr'q + h'^ 

= Sq-lU- WcTq - Aq (Pdiv(M • V)t + PdivgQ,(T, Vm))' 


Sq-lU- Vo-q - Aq {Vf^kdkdiVT)‘ 


- Aq{6im + diA ^dni)TQ kdkt 


Lm 


-Aq(dim + diA dm)diT'^^_^^,„M - A^(Pdivg„(T, Vm))', 

we bound Z^ez term by term. First of all, thanks to the commutate estimate Lemma 

10.25 in fli, we find that 


L 2 


S,_in • VMq - Aq (Pf„*5,diVT) 

qeZ 

(3.52) < C\\Vu\\l^ z z 2 ^(^i)||A,,divT||i 2 < Clln'^ll^ 


qeZ \q'-q\<4 
(d 


d 2 d 2 B 

^2,1 ^2,1 


2,1 


The second term contained in Z^ez^^'-z is nothing but ||(5™ + 5,-A ^d,n)fg^^^kdi.M'”\\ d_i, the 

estimate of which is easier. Indeed, using Proposition 12.11 and 12.21 directly yields 


Ms 


m„, + diA-^d^)Ts^,kdkT‘n\ 4-1 <C||Vn||i<»||VT|| d_,<C\\uW 

' r 2 n2 - 7 ^ 


J.m 


^ 2,1 


"^ 2,1 


"^ 2,1 


2.1 


Now we go to bound the third term. 


11(5™ + diA-%ddif'g^^..y \\. 4_, < Cll^^n'^ll.. < CWt'^^yW . d + CWt'^yW . d . 

^ 2,1 ® 2,1 ® 2,1 ® 2,1 

Using Proposition l2.2l once more and Bernstein’s inequality, it is easy to see that 
(3.53) liriyil . <C||5 ,t||b-, ||n'’|| 4 ,, < Clink’ll 4 ,J|t|| . . 

*^/t R 2 cx3,oo n 2 n 2 n 2 

•^ 2,1 ^ 2,1 '^ 2.1 '^ 2.1 

However, the estimate of \\T'^^^M\\ d is a little bit complicated. In fact, using Lemma [231 and Bern¬ 


stein’s inequality, we have 


^ 2,1 qeZ 


^ z z 2M2\\Aq{Sq^^2ddTAq,M)\y 

q-q'^\ 
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>112 


^ Z Z 

qeZ q-q'<l 

qeZ q-q'^1 q"<q' + l 

q^Z q-q'^l q"<q' + l 

< cZ Z Z 2('?"-^')5|2^"5 ||A,.t|L2 

qeZ q-q'<\ q"<q' + \ 


x(2?'(5 + l)||A./M"'||i2)" [29'5||A„,diV/,M' 


< CIIm^iZ ,||t|U . 

T + l 


4 , 

^21 


B 


2,1 


B 


2,1 


Finally, we would like to point out that the estimate of the term ||(PdivgQ,(T, Vm))'|| d _^ is somehow 


'^ 2,1 


“bad” due to the appearance of Y^\<k<d-\ \\Tri.i<dku‘^\\ d both for horizontal and vertical direction. 


Thus, there is no need to bound ||(PdivgQ,(T, Vm))'|| d_^ from different directions. In view of Proposi¬ 


tion |ZT] and Corollary 12.21 we easily get 


2,1 


||(PdivgQ,(T,VM)y|| rf 1 < C||gQ.(T, Vm)|| , 


(3.54) 


< CIIVull d\\T\\,d <C 




Hi 


• T 


Ki 


Combining the above estimates, we find that 
(3.55) 2]2^(^i)||/i'||^2 <C 

qsZ 

Substituting (13.511) and (13.551) into (13.481) yields 


+ llu^ll d _^ 


Hi 


• T"'" 

B 21 


"•7 


2.1 


rf 1 -F (1 - m)||cr'|| 1 -F (1 - m)||M'|| -Filer'll d^yd , 


2.1 


2,1 




< 


c 


(1 - m) 

/ 

f^l 

+ 


■j-i + rolL4-i 


". ^+1 + I|m II. 1+1 

Hi Hi 


HHi"" ■) 

i,i]M. idt' 


+ 


f 

f 

Jo 

C 




B}, B}, B}, B? 


2,1 


2,1 


2,1 


Hi 


di' 


^ 2,1 / 


, ft ||2 


2,1 


2,1 


[iIm'II 4 _i + ||cr' 

If .1 

'' 

dt' 

V Hi 

Hi ) 

J 


(1 - m )2 


o".^' + IrolU-i + 

^ 2,1 " 2,1 


I p 2 

2,1 ^ 2,1 


dt' 


^ 2,1 


/' 

■f 


rf_il|w^|| rf+, -F "-^"2 

Hi 


Hi 


5 !+* e^-' gj-' £3+* 

^ 2,1 ^ 2,1 ^ 2,1 ^ 2,1 y 

\ 


.4+1 + 
2.1 


till 2 


dul 


^ 2,1 


dt' 


d I dt' 


2,1 


for all / e { 1 , 2 , • • • ,d}. 
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Consequently, there exist two constants Ci and C 2 depending on d. Re and We, but independent of oj, 
such that (13.31) and (13.41) hold. 

The damping effect of r. 

In order to deal with the nonlinear coupling between u and r, it is necessary to take advantage of the 
damping effect of r. To do so, applying to the equation of t yields 

Tq 2a> ■ 

dtTq+Sq-lU ■VTq + —D{Uq) + Sq-lU ■ Vjq - Aq{u ■ Vt) - Aqga{T,Vu). 

Taking the inner product of the above equation with Tq, integrating by parts, we have 

2 a> 

= ^(D{Uq)\Tq) + (S q-lU ■ ^Tq - Aq{u ■ Vt)|t^) - {Aqga{T, Vm)|t^). 

It follows that 

\\Tq{t)\\i2 + :^l|TgllL;(L2) ^ C(||T^(0)||i2 + <j}2%Uq\\i2 + • Vt^ - Aq{t^kdkT)\\i2 


+ 11 

.d 


Multiplying the above equation by and summing over ^ € Z, we find that there exists a constant 
C 3 depending on d, Re and We, but independent of oj, such that 


IL + 


lUbI) 


< C3 || to ||.4 + a )\\ u \\ 


(3.56) < C 3 IItoII rf + m ||m' 


^2,1 


+ r 

'(bI^ ) Jo 


\\u\\,,J\T\\,,dt' 


LliBl ) 


+ M 


1 ■ T + l 


+ 


+\\u‘^\\,, )\\t\\ ,dt\ 

Jo b|i e|i J 


where we have used the following nonlinear estimates which can be obtained in a similar manner to 
those of (13.521) . (13.531) and (13.541) . 

Y 2^^IISq-lU ■ VTq - Aq(t^kd/cT)lli^2 < C||Vm||l<»||t|| rf , 

r J D 2 


qeZ 


and 


d < C||m|| rf_,il|T|| rf , l|gff(T,VM)|| rf < C||Vm|| rf IMI J J. 


2.1 


2,1 


2,1 


This completes the proof Proposition 13. II 


Remark 3.1. We can bound Zg€z 2^^2 without making a distinction between horizontal di¬ 

rection and vertical direction due to the appearance of the nonlinear term gafr, Vn). Our above 
estimates 0 /aim to reveal that the convection term u ■ Vt behaves better than 
ga{T, Vm). 

The rest part of this section is to close the estimates (I3.3I) - (I3.5I) in Proposition l3.ll Before proceeding 
any further, let us denote 

4 C 2 


Ad 

Aq .- 


aH ._ 
Aq .- 


fin := 


(1 - m)3 
I6C1C2C3 


(a'^(O) + 1 ), 

8 (C 2 + C3) 


(1 - 


■ exp 


(1 - 


(A^f (A''( 0 ) + mA''( 0 ) + R( 0 )), 


192CiC2c 2 n2(C^ + C 3 ) . , , V 

trrJr [ v-o,f ®®■ 
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Proposition 3.2. Assume that (u, r) is a solution of system (11.11 ) on [0, T\ satisfying 


(3.57) 

If 

(3.58) 

(3.59) 


A\T) < A[;, A\T) < A^, B{T) < Bq. 

(1 - 


4 


^A" + JA^^Ai + Bn < 


Bn 1 + 


0 0 u 2C2 
2 C 2 


Ag(A[; + Bo)< 1 , 


(1 - cof 


(^0 + So)) 


(l-w)3 


4Ci ’ 


< h 


(3.60) 

then there hold 

(3.61) 

Proof From (13.41) . it is easy to verify that 


2 C 2 C 3 


■A" < 


2“° “ ’ (1 -m)3 ° “ 2’ 


a\T)<'^-AI a\T)<^-A‘1, and B{T) <^-Bn. 


A\t) < 


C 2 


a\ 0) + A\t) (A\t) + B(0) + ^A'’(0l|n^lL^ . 


(1 - m)3 

+iAV)ll«^ll , +(aV) + 

2 ) 



LT(Blf) 


, +B(0l|n^|l 


(3.62) 
If 

(3.63) 


C 2 


(1 - m)3 


‘^(O) + A\0(a^(0 + B(0) + \^A\t) + ylA\t)Ad{t) + 5(0)A'='(t)|. 


4 


—An + -^IAqA^ + Bq < 


2 

(1 - 

2 C 2 ’ 


- (i-%3 ^ ^ ^ [0, r]. 


then it follows from (13.571) and (13.621) that 

(3.64) 

Furthermore, if 

(3.65) Ag(A{; + Bo)< 1, 
then (13.641) changes to be 

(3.66) 


A\t) < 


(a‘'( 0) + 1) - ^A^, te[0,r]. 


(1 -m)3 2 ‘ 

Next we go to bound B{t). To do this, combining (13.561) and (13.641) . we are led to 

5(0 < C3|5(0) + mA(S + mA^( 0 + f (||m'*|| 4 ,, + lln^ll 4 „)I|t|| ,dt'\ 

V Jo B|i B|i B|, ) 

< C 3 1 5(0) + mA* + 7 ^^^ (a"'(0) + (a'XO)^ + A\t)B{t) 


+ f(l|nl.4,, + 

Jo 


(1 - m)3 




2,1 

3 


"^ 2,1 


(3.67) < C,\B(0} + ju,A'i + 2^^. 


f 


A^(0)+ (Hull + 


■^ 2,1 


t|| rf JO I + -B{t), 

Hi 


provided 

(3.68) 


2 C 2 C 3 

(1 - ojf 


A'^ < — 
^0 - o' 
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Using Gronwall’s inequality, we obtain 

(3.69) 

where Bq is defined to be 

(3.70) Bo := 2Cs exp ( 2 C 3(1 + a;^)) (b( 0) + ^ 
To estimate we infer from (13.31 ) that 


Bit) < 2 C 3 exp ( 2 C 3(1 + A^)) (b( 0) + |mA[S + | < Bo, 

IC 2 C 0 

^OJ -T 

2 (1 - m)3 


A"'(0) 


A\t) < 


Cl 


(1 - m)3 




+ 


1 


1 M . £ + 1 


d ,dt' 


4Ci L/CB?, ) 2 Jo 




fill 


b|i 


d I dt 


2,1 


+||t|| .4 lli^^ll + 


L“(fi2j) ) ^0 ^211 


f 


llrll ddt' 


< 


Cl 


(1 - m)3 


a\ 0 ) + ||t|| . 


.. ,, ''+1 + 

Lr(Bl) LJbIj ) 


(3.71) 


+ 


73 . 


(1 2 


(1 - fa^)^ 
4Ci 
\ 


+ B(t) 


1 ■ T + l 

U (fi|i ) 


. ^+1 
^ 2.1 y 


d .dt'\. 

fill ' 


In view of (13.641) . we have 

,d 


2 C 2 

U II rf < -T 

(1 - w)J 


(a‘'(0) + A^(0 (a\0 + Bit))) Bit) 


(3.72) 


< . A^(0)g(t) + 77 ^^^ (a\ 0 + B(t)) Bit)A\t). 


(1 - m)^ (1 - m)^ 

Substituting (13.721) into (13.711) . and using (13.691 ). one deduces that 


A\t) < 


U\0) + . a\0)Bo + (Bit) ( 1 + 


(1 - m)3 

"73 


+ 


(3.73) 

If 

(3.74) 

then (13.731 ) reduces to 


2 B 


.i+i + 
2 

2,1 


(I - m)3 


Cl 


(I - m)3 


(A\t) + B(t))JJ 


A ^(0 


rA“ + ^ 


Bn I + 


(I -m)3“° ■ 2 
2 C 2 

(1 - 01)3 


fi£i 


2.1 


d-idt + —A^*(?). 




(l-cu)^ 

4Ci ’ 


a7o < 


2Ci 


(3.75) 


(1 - m)3 
"73 

+ 


A'JO) + 


2 C 2 


(1 - m)3 


A70)Bo 


Cl 


j I + I- tA^ + — 

B|;' \(l-m)3 « 2 


Gronwalll’s inequality implies that 


(3.76) A7t) < 
provided 


2Ci 


exp 


(l-m)J "\(l-m) 


2Ci 


K1+A71 + 


Cl 


(1 - m)3 


dj\\u^\\d_,dt' 


2 C 2 


A^)) A70) + 


(1 - m): 


rA70)Bo , 


-A^ < 1 

2 0 - 


(3.77) 
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Substituting (13.701) into (13.761) . we find that 

2C^ 


A'\t) < 


2Ci 


(1 - cu)3 


exp 


i^(l+A^)2j|A''(0) 


4 C 2 C 3 

+ —--7 exp 


(3.78) 


< 


< 


(1 - ujf 

2Ci 


{l-OJ, 

(2C3(1 + A^)) (a^(O)B(O) + a;A^(0)A9 


(1 - a)f 


exp 


8 C? 


8C1C2C3.( 

^8(C2 + C 3 ) 

(1 -m)6 

. (l-m)6 


4 C 2 C 3 
1 - W)3 

(A;!)d(A'‘(0) + (ja‘'( 0) + B(0)) = f aJ. 




Finally, substituting (13.781) into (13.691) . we get the estimate of B{t) in term of A*(0) + a;A'^(0) + B(0): 

^ 8 (C 2 + C3). . 2 ) 


48 C 1 C 2 C, 

®(0 ^ - -exp 


(3.79) 


(1 - m )6 

“* 2 ^^ 
.2^3 


(4C3Ag 


exp 


(1-m) 


^ (aQ j (2B(0) + a'*( 0) + 2mA^(0)) 


96CiC^C^ 


(1 - m)^ 


■ exp 


12 (C^ + C3) 

(l-m )6 


(aQ'] (a^O) + mA''(0) + B(0)) = jBo- 


The proof of Proposition I3.2l is completed. 


4. Proof of Theoreiv[TA] 

Now we are in a position to prove the maim result of this paper. 


Proof. The local existence and uniqueness of the solution (u, r) € £|. to system (11.11) with initial data 

B^i ^(^21 been proved by Chemin and Masmoudi fTl. Let T* be the lifespan 

of (m, t) obtained in 0- Define Ti be the supremum of all time T' € [0, T*) such that 

(4.1) A''(t) < A^, A‘^{t) < A^, B{t) < Bo for all t e [0, T']. 

Choosing Cq satisfying (11.81 ) so large that (I3.58I )- (I3.60I ) hold, then thanks to Proposition 13.21 for all 
0 < t < Ti, we have 


(4.2) aV)<^A|S, A^(0<^A^, B(0< ^Bo- 

According to standard continuation method, we obtain Ti = T*. Moreover, (14.21) and the embedding 
B^i ^ L°° imply that 

(4.3) r (||VM(t, •)||ioo + ||T(f, •)||z,~) dt < 00 , 

Jo 

provided (11.81) holds. Thus, the blow-up criterion in 0 enables us to conclude that T* = 00 . This 
completes the proof of Theorem ll.il □ 
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